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Equilibrium free energy measurement of a confined electron driven out of equilibrium 
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We study out-of equilibrium properties of a quantum dot in a GaAs/AlGaAs two-dimensional 
electron gas. By means of single electron counting experiments, we measure the distribution of 
work and dissipated heat of the driven quantum dot and relate these quantities to the equilibrium 
free energy change, as it has been proposed by G. Jarzynski [Phys. Rev. Lett. 78, 2690 (1997)]. 
We discuss the influence of the degeneracy of the quantized energy state on the free energy change 
as well as its relation to the tunnel rates between the dot and the reservoir. 


Equilibrium thermodynamics is a fundamental branch 
of physics providing tools to make predictions of macro¬ 
scopic many-particle systems independent of detailed mi¬ 
croscopic processes governing their properties. In the 
recent trend towards smaller systems, which deviate 
strongly from the thermodynamic limit, fluctuations de¬ 
parting from the equilibrium state often become promi¬ 
nent and non-equilibrium dynamics needs to be taken 
into account. The discovery of fluctuation relations mm 
and their experimental tests |3]-[T2] are major steps to¬ 
wards understanding the evolution of small systems down 
to the atomic level. Elementary building blocks for ther¬ 
modynamics on a microscopic level are discrete quantum 
states which are, for example, used for the logical ele¬ 
ments of quantum bits. 

Here, we study a single discrete energy level in a quan¬ 
tum dot coupled to a single thermal and electron reser¬ 
voir. Driving the quantum dot out of equilibrium with 
respect to the reservoir allows us to demonstrate exper¬ 
imentally the connection of an equilibrium quantity, the 
free energy, to the non-equilibrium dynamics as predicted 
by the Jarzynski equality [T]. The Jarzynski equality has 
been tested in systems involving many energy levels and 
for the special case of cyclic drive protocls [4]. Other 
experiments exist, where the Jarzynski equality has been 
used to determine the equilibrium free energy change in 
systems where precise calculations thereof are difficult 
to obtain ISIS]. In this experiment, we drive a single 
electron occupying a fully quantized energy level and we 
extract the equilibrium free energy along the full drive 
trajectory. The results are consistent with the theoret¬ 
ical predictions from equilibrium thermodynamics. We 
analyze the influence of a degeneracy of the discrete en¬ 
ergy level on the free energy. We further discuss gen¬ 
eral limitations of fluctuation relation based experiments, 
and particularly the importance of quantum dot (QD) 
systems as testbeds for statistical physics on the single¬ 
particle level. 

The ability to measure the free energy change in any 
process allows estimating the probability for its sponta¬ 
neous occurrence. Even more, the free energy is a direct 
measure of the maximum work that may be extracted 
from a process or the minimum work which needs to be 


invested to activate it m- A good characterization of a 
system in terms of the free energy is furthermore impor¬ 
tant for experiments where drive is applied in order to 
cool the system, or to retrieve information about system 
properties, as well as for Szilard’s engine and Maxwell’s 
demon type experiments. Our results for the elementary 
single level system paves the way for understanding small 
systems with a more complex state composition, which 
is the case, for example, in molecular structures coupled 
to a thermal bath laiHiis]. 



FIG. 1. (a) Scanning electron micrograph of the device 

formed on a GaAs / AlGaAs heterostructure with a 2DEG 90 
nm below the surface (dark grey area). Ti / An gates (bright 
fingers) are patterned with e-beam lithography to capacitively 
deplete the 2DEG. The QD is formed in the place indicated 
by the blue circle and is tunnel coupled to a reservoir, as 
shown by the red arrow, while the other barrier is fully closed. 
An arbitrary waveform generator (AWG) is connected to the 
plunger gate to drive the QD. (b) An energy diagram indicates 
a tunneling process between the QD and the lead reservoir at 
energy E measured from the Fermi energy /i. 

As a first experimental step we characterize the QD- 
reservoir system in thermodynamic equilibrium. A typ¬ 
ical time trace of the charge detector current is shown 
in Fig. [^a). In our experiment, we form a QD in a 
GaAs / AlGaAs two-dimensional electron gas (2DEG), 
as shown in Fig. The confinement can be controlled 
by applying negative voltages to top-gates. The voltage 
VpG applied to the plunger gate tunes the energy level 
E linearly. The QD is coupled to a large contact re¬ 
gion of the two-dimensional electron gas acting as the 
electron and heat reservoir characterized by the temper- 
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ature T ^ 40 mK. We measure the tunnel rates for the 
last electron uniiTT] with single-electron counting tech¬ 
niques, utilizing a nearby quantum point contact as a 
charge detector [Ulllg] and find electron tunneling at 
frequencies of the order of 100 Hz. 

The occupation probability / = 4n/J2f=ii4n + 

^out) and tunneling rates Win/out = N/ i-n/out’ 
the number of tunneling events N, are estimated from 
similar time traces taken at the indicated plunger gate 
energies p^H2Q] and presented in Fig. [^b). The energy 
relaxation time for electrons injected from the QD into 
the Fermi sea is assumed to be much shorter than the 
time scales relevant for electron tunneling between the 
QD and the reservoir, and the chosen drive frequencies 
of less than 10 Hz [21]. Electron tunneling between the 
QD and the lead is described by Fermi’s golden rule, 
with energy independent tunnel coupling constants Fin 
and Tout fhe Fermi distribution f{E,T) describing 
the lead occupation at an energy E measured from the 
Fermi energy of the lead, fi. The resulting tunneling rates 
are 

fHin(^) = rin/(F;,T) (la) 

f^out(^) = rout[i-/(^,T)]. (lb) 

From a weighted least-mean-square fit of the measured 
tunnelling rates to Eqs. 0 . the tunnel coupling con¬ 
stants, the position of the electrochemical potential, ji = 
0, used as the zero-energy reference, as well as the tem¬ 
perature in units of plunger gate voltage of the lead 
are extracted. The fit is shown as solid blue line in 
Eig. I^b). In the same figure, the solid green line is a 
least-mean-square fit of the estimated occupation prob¬ 
ability / to a second Eermi function, with the tempera¬ 
ture and the resonance as fit parameters. Erom this fit, 
we recover the same temperature. However, we find that 
the occupation probability is one half around an energy 
value which is offset with respect to the Eermi energy by 

E/kT = om±om. 

Erom the partition function of the QD-reservoir sys¬ 
tem one can see that the offset of the energy where 
/ = 1/2 with respect to the Eermi energy is determined 
by E = /cTln(d) in the case where one electron is filled 
into an empty d-fold degenerate energy state. Eor the 
last electron, d = 2 due to spin and hence, the occupa¬ 
tion probability is one half at E = /cTln(2), as found in 
the experiment. At E = /i, the probability for the energy 
level to be occupied is twice the probability for the state 
to be unoccupied. As shown in the inset of Eig. [^b), 
the system obeys the detailed balance condition for the 
tunneling rates, bFin/kFout = dexp(—E//cT), with d = 2, 
as found from the measured d = Fin/Tout = 1.97 ± 0.04. 

After having characterized the QD-lead system in equi¬ 
librium, our goal is to measure the equilibrium free en¬ 
ergy change AE of the QD between an initial and a final 
state differing by a value AE oc —Vpg while it is driven 



FIG. 2. (a) A typical time trace of the charge detector 

current (Icd). (b) Tunnel rates (filled blue for tunneling 
in and open blue dots for tunneling out) and relative fre¬ 
quency of occupation (green squares, filled for occupied and 
empty for non-occupied) extracted from CD time traces of 
60 s. The abscissa is the negative plunger gate voltage mea¬ 
sured from the Fermi energy and scaled by the temperature, 
T = 223 /iV in units of Fpg- The error bars indicate statisti¬ 
cal errors assuming Gaussian fluctiations. The solid blue lines 
are weighted least-mean-square fits to Eq. Q, from which 
Fin = 41.8 ± 1.6 Hz and Tout = 21.2 ± 0.6 Hz and the temper¬ 
ature are determined. The dashed green lines are fits to Fermi 
functions centered around AE = ksT h\{2). Inset: The ratio 
of the two tunnel rates (black circles) is plotted as a function 
of energy from /x, together with the exponential behaviour as 
expected from the detailed balance condition (red line). 


out of equilibrium with respect to the reservoir by ap¬ 
plying a time-dependent plunger gate voltage Epc- The 
Jarzynski equality [l] provides the necessary tools for the 
theoretical description of such an experiment. It states 
that for an arbitrary drive protocol, the performed work 
AW during the drive satisfies 



where the average is taken over different repetitions of 
the drive protocol and AE is the equilibrium free energy 
change between the initial and final state of the drive. As 
the work can be evaluated for non-equilibrium situations, 
the Jarzynski equality allows to determine the equilib¬ 
rium quantity AE with a non-equilibrium measurement. 

In our experiment, we apply a sinusoidal voltage ramp 
kpG(Q to the plunger gate, which changes the energy of 
the electronic state in the QD from Eq to Ei and per¬ 
forms work on it as long as it is occupied [see Fig. [Tib)]. 
The drive protocol, as shown in the inset of Fig.^ is 
characterized by the frequency / and amplitude A of 
the sinusoidal part of the protocol. Together, A and / 
determine the maximum steepness of the ramp and the 
initial (Eq = —A) and final (Ei) values of the QD en- 
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FIG. 3. (a) Examples of drive realizations in upward (blue) 

and downward (green) direction, plotted as solid lines. The 
applied drive signals are shown as dashed lines. The driving 
is done by applying a voltage Vpg oc — E to the plunger gate 
through an arbitrary waveform generator. From the position 
of the tunneling events (grey dashed lines) the heat and work 
of each realization are extracted. Bottom: The distribution 
for AQ (left) and AW (right) measured from approximately 
20000 realizations. Blue circles originate from upward and 
green squares from downward drive. The solid lines are rate 
equation calculations without free fitting parameters. 

ergy level. We distinguish two drive directions, upward 
{El > Eq) and downward {Eq > Ei), and employ a 
waiting time of At = 0.5 s for equlibration in-between. 
The respective initial and final states, g^j^d 

^iUp,Down) ^ are indicated in the figure. Typical time- 
traces showing the charge-detector signal for four indi¬ 
vidual voltage ramps are shown in Fig. |^a). In most 
cases, as shown on panel I, when the energy level of the 
QD starts well below the Fermi level and is raised above 
the Fermi level, the electron leaves the QD. If the drive 
frequency is large compared to the tunneling rates, re¬ 
alizations without tunneling events become probable, as 
shown in panel II. For each realization, we determine the 
work AW perfomed on the QD. When driving an occu¬ 
pied energy level of the QD, work is performed on the 
electron by the voltage source amounting to 



where n{E) G {0,1} denotes the occupancy of the QD, 
which is directly measured by the charge detector sig¬ 
nal. For example, in realization I, with the drive ap¬ 
plied symmetrically around zero and A = 0.45 kT^ the 
work performed on the electron in the QD is AW = 
(0.45 + 0.11) kT = 0.56 kT. Note that the maximum 
amount of work that can be performed on or by the QD 
system is given by ±2^4. 

In addition, for each realization we determine the heat 
AQ dissipated in the contact from the energies Ei of the 


QD state at which tunneling events occur. Each tunnel¬ 
ing process contributes AQ^ = Ei for tunneling out and 
A^Qi = —Ei for tunneling into the QD. For example, in 
realization I, the electron tunnels at 0.11 kT above /i, 
hence Q = 0.11 kT due to the relaxation of the electron 
in the lead. If more than one tunneling event occurs dur¬ 
ing a single realization, as shown in panels III and IV, 
the individual processes are added up, 

AQ = Y,AQi = Y,EiSi, (4) 

i i 

where Si distinguishes between tunneling-out processes, 

= 1, and tunneling-in processes, Si = —1. The maxi¬ 
mum amount of heat which can be dissipated in the lead 
during a single drive realization amounts to ±2^4. 

With the analysis described above, we now determine 
the relative frequency distributions of AW and AQ as 
shown in Fig.j^b). The distributions for driving upwards 
are plotted as blue circles, while green squares are used 
for the distributions resulting from driving downwards. 
With the tunnel rates given in Fig. |^b). Win = 42 Hz 
and Wont = 21 Hz, and a drive frequency of / = 10 Hz, 
realizations as shown in panel H of Fig. [^b), i.e. with¬ 
out tunnel events and hence zero dissipated heat, are very 
probable, leading to a prominent delta-peak at AQ = 0 
in Fig.j^b). In these realizations, no work is done on the 
QD if the level stays empty during driving, leading to a 
peak at AW = 0 in Fig. [^b). If it is occupied, the work 
equals plus or minus twice the drive amplitude for up¬ 
ward or downward drive direction, respectively. Hence, 
we observe two additional peaks, at AW = ±2H, also 
corresponding to the single delta-peak at zero dissipa¬ 
tion. From Fig. [^b), it is apparent that the relative 
frequency for |Q| < H is much higher than for |Q| > A. 
This is due to the fact that the maximum dissipation in 
a single realization depends on the occupancy of the ini¬ 
tial and final state. For example, when driving an energy 
level downwards, the dissipated heat can only exceed A 
if the state has been occupied in the beginning, as in 
realization IV. Contrary, AQ < H is only possible for 
downdrives if the final state is the less probable unoc¬ 
cupied state. Since the QD is thermally equilibrated by 
the waiting time before each drive, the occupancy of the 
initial state is given by the Fermi distribution and re¬ 
alizations with large |Q| are suppressed. This effect, to¬ 
gether with a comparably fast sinusoidal drive, also leads 
to the peaks found at Q = ±A: they are mostly due to 
realizations starting in the more probable initial state 
and having a single tunnel event at the very end of the 
drive. Generally, the sharp features found in the proba¬ 
bility distributions shown here are the signatures of non¬ 
equilibrium dynamics. For slower drives, the probabil¬ 
ity for intermediate dissipation |Q| ^ 0 increases and we 
find more Gaussian shaped distributions centered around 
zero dissipated heat (data not shown). This corresponds 
to the adiabatic limit, where a Gaussian shape with a 
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width given by the fluctuation-dissipation theorem is ex¬ 
pected [22] , as discussed in more detail in Ref. [4] and [10] 
The detailed shape of the distributions depends on the 
relative amplitude A/kT and frequency f/Tin^out of fho 
sinusoidal drive. For example, the distributions for work 
and dissipated heat are asymmetric with respect to the 
two drive directions because of the different rates for tun¬ 
neling in and out of the QD, as discussed above. Using 
the standard rate equation approach |4] with parameter 
values extracted from the measurements of Fig. [^b), we 
calculate the probability distributions for the work and 
the dissipated heat in each drive direction and find very 
good agreement with the experimental data. 



FIG. 4. The free energy of the QD during the drive, obtained 
with non-equilibrium measurement, measured with respect to 
the Fermi energy. Experimental results are plotted in green 
for the downward and in blue for the upward drive direction. 
Different markers denote independent experiments with dif¬ 
ferent drive parameters. The theory curve, Eq. is plotted 
as solid red line, for d = 2. Dashed grey lines indicate the 
expected linear increase at low Ei as well as the saturation 
value for Ei > Ep- Inset: One full period of the drive is 
plotted in blue. A waiting time of 0.5 s ensures thermal equi¬ 
libration of the initial state. Green dashed lines indicate the 
values of Ei for which we extract AF from the experiment. 

Utilizing Eq. we now determine the equilibrium 
free energy change of the QD from the statistics of per¬ 
formed work. The prerequisite that the initial state must 
be thermally equilibrated [1] is fulfilled by the waiting 
time. We determine the free energy change between the 
fixed initial and any freely chosen final state within the 
drive by analyzing AW as described above. As a result 
we get the free energy difference of the initial state and 
the chosen state within the drive. By varying the final 
state in the analysis, the change in free energy is deter¬ 
mined along the full drive trajectory, as shown in Fig.|^ 
Different markers are used for the three measurements 
performed with the same gate configuration and QD res¬ 
onance but with different values for A and /, as indicated 
in the figure. Again, blue color denotes upward and green 
denotes downward drive direction. Each of the six drive 


protocols gives an independent estimate for AF. From 
the extracted AW we also find that the average work (W) 
always exceeds the free energy change, as is expected for 
a system which is driven out of equilibrium. 

Theoretically, the free energy is calculated from the 
difference of the partition sums of the initial (i = 0) and 
final (i = 1) states [13], 

AF = kTln(Ay Zi = l+de^p(^^'^ (5) 

The theoretical curve is plotted as a solid red line in Fig.[^ 
for d = 2, which is the value obtained from the equlib- 
rium characterization above. Compatible with theory, 
we observe a linear increase of the free energy in the first 
part of the drive, where the electron is lifted to an energy 
level still well below the Fermi energy. The increase in 
AF slows down around the Fermi energy, as the proba¬ 
bility of the QD to be occupied decreases. At energies 
well above the Fermi level, the occupation probability 
approaches zero and the free energy saturates. The six 
independent measurements of AF agree well with the 
theoretical value along the whole drive trajectory. 

Small systematic deviations from the theoretical value 
of AF are found only for the data resulting from up drives 
at / = 2 Hz and A = 4.03 kT. Experimental limitations 
are given by the stability of the sample as well as the 
bandwidth of the set up. The bandwidth limits the ac¬ 
curacy in the determination of the exact tunnel time in 
drive schemes where large amplitude and frequency are 
combined. These errors lead to uncertainties in AW and 
thereby also in AF. The sample stability, on the other 
hand, sets a limitation to the number of repetitions which 
can be performed in a given configuration of the QD. 
Most importantly, the drive must be applied around the 
chosen operation point, equally for every repetition. En¬ 
ergy drifts in time therefore lead to systematic errors. 
Although we observe stable sample configurations over 
many hours, we implemented a feedback mechanism in 
order to align F = 0 to the Fermi level /i every quarter 
hour. Offsets are taken into account in the analysis. 

As demonstrated in Fig. we have found good agree¬ 
ment between the experimental determination and theo¬ 
retical expectation of the free energy change in our well- 
characterized system. This supports the idea that the 
free energy can be measured at any point of the drive 
in a system far out of equilibrium, independent of the 
drive parameters, by utilizing the Jarzynski equality. In 
the original work by C. Jarzinsky, an experimental test of 
equation!^ is suggested in a nanoscale system weakly cou¬ 
pled to a reservoir, where the fluctuation of the work is 
less then kT and which evolves deterministically under its 
Hamiltonian [T] . Even though the phase space evolution 
of the QD in our case is intrinsically non-deterministic 
due to tunneling events and the fluctuations are of the 
order of kT drive amplitudes, we find Eq. to be valid 
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and the QD system to be a very suited testbed for the 
Jarzinsky equality. 

As known from equilibrium thermodynamics, the free 
energy is a state variable and therefore independent of the 
trajectory. When the system is driven out of equilibrium, 
such trajectory-independent quantities are not generally 
available and the situation becomes more complex. In 
this work we have shown that an equilibrium quantity, 
the free energy of a two-fold degenerate single electron 
state, can be measured with an experiment driving the 
system far from equilibrium. Our results show that ther¬ 
modynamic quantities such as work and dissipation can 
be understood on the level of individual quantum states 
and in nonequilibrium. The good agreement with the¬ 
oretical calculations proves that the simple system we 
have used suits well for studying more complex phenom¬ 
ena out of equilibrium. As a direct consequence of the 
presented work, we expect the relation between thermo¬ 
dynamics and information theory as well as influences of 
dissipation on future fast drives of qubits to be studied 
in quantum dot systems. 
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